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^ ' ABSTRACT: The complementarity between dark matter searches at colliders and in under- 



ground laboratories is an extraordinarily powerful tool in the quest for dark matter. In 



the vast majority of the analyses conducted so far these dark matter detection strategies 
r S^ . have been profitably combined either to perform global fits in the context of certain particle 

physics models (e.g. the CMSSM) or to estimate the prospects for a direct dark matter 
detection given the LHC potential of discovering new physics beyond the Standard Model. 
In this paper we propose an alternative strategy to combine direct and collider dark matter 
searches: employing the potential of the upcoming generation of 1-ton direct detection ex- 
periments, we show that for certain supersymmetric configurations it is possible to translate 
the information encoded in an hypothetically discovered direct detection signal into classes 
of expected signals at the LHC. As an illustrative application of our method, we show that 
for a 60 GeV neutralino thermally produced via resonant annihilations and identified by 
a 1-ton direct detection experiment, our approach allows to forecast a clearly identifiable 
prediction for a LHC final state involving three leptons and missing energy. The strategy 
presented in this paper to systematically translate a direct detection signal into a predic- 
tion for the LHC has the potential to significantly strengthen the complementarity between 
these two dark matter detection strategies. 
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1 Introduction 

The nature of the dark matter (DM) component of the Universe remains unknown. One 
of the most popular solutions to this puzzle is the scenario in which DM is made of a 
beyond-Standard-Model (SM) particle which is stable, in thermal equilibrium in the early 
Universe and non-relativistic at chemical decoupling; the relic density for such a state 
scales approximately with the inverse of its pair annihilation rate into lighter SM particles 
and matches the abundance of DM as determined from cosmological measurements for 
a weak-interaction-strength coupling: this is the essence of the celebrated paradigm for 
weakly interacting massive particles (WIMPs) as DM candidates. On top of the fact that 
there are numerous extensions of the SM naturally embedding a WIMP DM state, to some 
extent the popularity of this framework relies on its very rich phenomenology, with a set of 
complementary detection strategies stemming from the framework itself on rather model- 
independent bases: One in general expects a tight correlation between the annihilation rate 
of DM in the early Universe (at a temperature of about one twentieth of the particle mass) 
and that in dark halos today (at relative velocities as small as 10 to 1000 km s _1 ) giving 
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rise to a non-standard astrophysical source of photons, leptons or antimatter through which 
indirectly detect the DM particles. On the other hand, applying CP symmetry arguments 
to the annihilation process, one expects to produce DM states in inelastic scatterings of SM 
particles at high energy colliders, such as the LHC. Finally, based on crossing symmetry, 
one generally predicts that the weak coupling to ordinary matter in the early Universe 
corresponds to a small but still significant elastic scattering cross section for DM particles 
on ordinary matter, through which directly detect the DM WIMPs forming the halo of the 
Milky Way. Direct, indirect and collider searches for WIMPs have then been primary goals 
in the experimental effort towards the identification of the nature of DM in the latest twenty 
years or so, having reached today the stage when, in all three techniques, experiments have 
the capability to single out or exclude viable DM models. 

In this paper we investigate some aspects of the complementarity between direct de- 
tection (DD) and collider searches. The goal is to investigate the potential of the upcoming 
generation of DD experiments, taking as reference a setup resembling the 1-ton realization 
of the XENON detector, and cross correlate the information encoded in a hypothetically 
discovered signal with the expected signals at the LHC. Modeling the coupling of DM and 
SM states at an effective operator level, such as in the recent analysis of, e.g., [1] one can 
drive a direct link between the two detection strategies; on the other hand, when consid- 
ering explicit WIMP models cross symmetry arguments rarely apply, with the couplings 
ensuring the appropriate relic density most often not corresponding to those setting the 
WIMP-nucleus elastic scattering cross sections and those responsible for the production of 
DM states in LHC's pp collisions. We choose then here to focus on a definite scenario, 
assuming that the DM particle is the lightest neutralino in the minimal supersymmetric 
extension to the SM (MSSM). This choice is motivated by the fact that, with a number of 
degrees of freedom which is more than twice those of the SM and the large number of free 
parameters involved (we will refer to a simplified version of the so-called phenomenologi- 
cal MSSM, or pMSSM, allowing for 12 free parameters), this framework is general enough 
to embed several different mechanisms to provide the correct relic abundance and/or a 
scattering rate at a detectable level. At the same time, the neutralino is the most closely 
investigated WIMP DM candidate within the most widely studied framework for searches 
for physics beyond the SM at colliders (as well as in rare flavor processes): we will then 
exploit the several well-known patterns for this DM state as criteria to explore the com- 
plex parameter space and interpret results, as well refer to a set of well-established search 
channels to extract guidelines which are actually more general than for the MSSM itself. 

Our analysis starts with the choice within the MSSM of a set of few benchmark mod- 
els representative of WIMP DM candidates which, in the plane spin-independent elastic 
scattering cross section per nucleon versus DM mass, appear below the region excluded by 
DD searches - with the most relevant limit stemming from the null result recently reported 
by the XENON100 experiment [2] - but still giving a sizable number of scattering events 
for a 1-ton detector. All benchmarks are selected checking that the lightest Higgs boson 
in each model matches the Higgs- like state recently discovered at the LHC, and fulfilling 
current constraints on beyond-SM states from rare process like b — > sj and B s — > as 
well as general trends from direct SUSY searches at colliders. For each benchmark, first a 
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set of mock data for the recoil spectrum in a 1-ton detector is simulated; such spectrum 
depends, on top of the spin-independent scattering cross section and the DM mass, on the 
phase space distribution of DM particles. Uncertainties on the latter have been recently 
investigated in Refs. [3, 4], implementing a large set of dynamical observables for the Milky 
Way and under the simplifying (but reasonable) approximations of spherical symmetry and 
isotropy; we will exploit results of such analysis in our approach here (complementary meth- 
ods to study the local DM distribution have been discussed for instance in Refs. [5-10]). 
Our simulated DD data samples consist in a number of measured recoil energies varying 
from approximately 30 to 190, depending from the benchmark model under analysis. Before 
analyzing these data, a statistical test has been performed on the simulated data samples 
to assess whether they provide a faithful realization of their true underlying distribution. 
This aspect has been addressed with great care since statistical fluctuations in the distri- 
bution of the observed recoil energies can potentially affect the DM mass and scattering 
cross section determinations significantly, in particular for heavy DM candidates with low 
scattering cross sections [11]. Mock data are then analyzed employing a Bayesian approach 
supported by a Markov Chain Monte Carlo (MCMC) Likelihood scanning technique of the 
parameter space defined by both the particle physics model parameters and the parameters 
related to the DM phase space distribution function. The results of this analysis regarding 
the preferred values for the DM mass and scattering cross section will be presented in terms 
of Bayesian credible regions. Having derived through this procedure the accuracy within 
which the relevant DM properties, namely mass and scattering cross section, are recon- 
structed and their mapping on the MSSM parameter space, we take this input to outline 
the most suitable collider search strategies and argue next future detection prospects. This 
second part of our analysis has been performed by focusing on the regions of the DM para- 
meter space characterized by the highest statistical weight (i.e. with the largest Bayesian 
posterior probability density resulting from our DD analysis) and simulating for a subset 
of models, sampled within these regions, the ATLAS detector response in its configuration 
with 7 TeV center of mass energy. This allowed us, in particular in the case of one of the 
benchmark investigated here, to effectively translate the simulated DD signals studied in 
the first part of this work into definite predictions for certain LHC observables. This proce- 
dure defines a method to efficiently exploit the upcoming generations of 1-ton DD detectors 
to forecast valuable predictions for the LHC and it constitutes the main achievement of 
this paper. The paper is organized as follows: in section 2 we provide a description of the 
benchmarks studied in this paper and review the criteria which brought to their selection. 
Section 3 is instead devoted to a study of the capability of a ton-scale DD experiment in 
reconstructing mass and spin-independent cross section of the DM candidates associated 
with these benchmarks. In section 4 we will translate these informations on the DM mass 
and scattering cross section into a prediction for certain classes of LHC final states while 
our conclusions will be summarized in section 5. 
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2 Determination of the benchmark points for Direct Detection 

This section is devoted to the determination of the benchmark points for the DD analysis. 
We will start by introducing the MSSM considered in this work and then we will review 
the dependence of the scattering cross section from the SUSY parameters stressing the 
implications in our framework. 

2.1 pMSSM setup and general constraints 

For our investigation of the MSSM parameter space we will pursue a phenomenological 
approach adopting the framework dubbed pMSSM [12]. Its most general realization features 
22 parameters defined at the Electro- Weak (EW) scale: 

• Mi, M.2 and M3 parameters corresponding to the gaugino masses; 

• \x, rriA and tan/3 parameters describing the Higgs/Higgsino sector; 

• sfermion soft-SUSY breaking mass parameters rriq, rriu R , , mf , m^, equal for 
first and second generation ; 

• trilinear soft terms A u , A& and A e equal for the first and second generation; 

• third generation soft mass parameters rriq 3L , fni R , mi R , m f 3i ) m f R 'j 

• third generation trilinear soft terms At, A\> and A T \ 

with the further assumption that no additional sources, with respect to the SM, of flavor 
and CP violation are present. The exploration of such large parameter space is, however, 
computationally very demanding; for this reason we have considered a few limits reducing 
its dimensionality, paying attention however to still cover all the scenarios that are, possibly, 
phenomenologically relevant. 

First of all we have neglected contribution to all observables from squarks of the first 
two families fixing their mass at about a few TeV, decoupled from the rest of the spectrum. 
In addition, their corresponding trilinear A— terms have been set to zero. This choice is 
motivated by the fact that current null results of SUSY searches at LHC constraint more 
severely masses of the squarks of the first two generations, as well as the gluino mass, ruling 
out values even about 1 TeV [13, 14] 12 . For the same reason we have considered only values 
of the gluino mass greater than 600 GeV. 

For simplicity we have also further reduced the parameter space in the sfermion sector 
by assuming mj = = mg fl , = Trif 3L = m? R , and rriq. iR = = mi R . We have 
finally also set A^ = A T . 

1 Most of the exclusion limits available in the literature rely on oversimplified models, mostly on the 
CMSSM, while their translation into a multi-parameter model like the pMSSM is not straightforward. In 
the scans considered in the next subsection we will not apply excessively tight constraint determining more 
accurate limits case by case through our detector simulations. 

2 Customarily considered collider searches can be evaded in the case of spectra enforcing particularly 
low mass splitting among particles involved in the production as subsequent decay chains, thus weakening 
exclusion limits even for first family squarks. We will neglect however this possibility in order to avoid 
dealing with a parameters space with an excessively high dimensionality. 
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Given this setup we have performed a series of flat scans in order to determine the 
relevant scenarios for the Monte Carlo study. The 12 remaining SUSY parameters have 
been varied over the following ranges: 

10 < Mi < 3000 GeV 

100 < M 2 ,n < 3000 GeV 

600 < M 3 < 3000 GeV 

100 < mg 3L , mg 3R , mj, < 3000 GeV 

90 < m A < 3000 GeV 

3 < tan p < 60 

-9 < A t ,A ft < 9TeV (2.1) 

SUSY spectra have been generated with the package SUSPECT [15] and required to be 
physical, i.e to satisfy the EW breaking condition, and to have a neutralino as LSP. We 
have in addition considered a series of general bounds. 

We have first of all required that the generated spectra provide a CP-even Higgs boson 
mass in agreement with the recent LHC discovery [16] of a neutral boson of mass of around 
125.5 GeV and couplings compatible with the ones of the SM Higgs boson (there are slight 
indications, not statistically decisive, of some deviations from the standard values, see for 
details [17, 18]). The MSSM setup most likely in agreement with this result is the so-called 
decoupling regime [19] in which the light CP-even Higgs state h has SM-like couplings and 
is required to have mass in agreement with the experimental determination, while all the 
other Higgs states are degenerate in mass and lie at a sensitively higher scale. In this regime 
the Higgs mass m,h is approximatively given by the sum of tree-level term and stop-top loop 
contributions: 

"* " mz cos 213 + ^ ( log (f ) + § (' - iw) ) (2 - 2) 

where Xt = At — /xcot/3, namely the stop mixing parameter, and M s = *Jrn% 1 with 
t\ and t 2 being the two stop mass eigenstates. Matching the experimental value of the 
Higgs mass is rather challenging for the MSSM. Indeed it requires a large contribution from 
loop-corrections which achieve the correct value only in the case the stop mass scale lies 
above around 3 TeV unless the stop mixing parameter Xt = At — fi tan /3 fulfills the so called 
maximal mixing condition [20, 21], namely Xt ~ y/6M s 3 . In this work we have imposed 
the requirement 123 < < 129 GeV, taking into account the experimental uncertainties 
as well as the ones affecting the numerical computation. 

We just mention that correct values of the Higgs mass can be also achieved by config- 
urations apart from the decoupling regime, featuring light, i.e. close in mass to the LHC 
candidate, masses for the other Higgs states (see e.g. [23] for a review). However these are 

3 Milder constraints are obtained in the minimal extension of the MSSM, the so called NMSSM, in which 
the Higgs couples with and additional superfield which provides an additional one-loop contribution to its 
mass [22]. 
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severely constrained by the searches of CP-odd Higgs decays into tau pairs [24, 25] and, as 
we will clarify in the following, by the limits from Br{B s — > /x + /z~). We finally notice that 
light DM scenarios may influence Higgs phenomenology allowing for the decay h — > XX- 

The invisible decay width of the Higgs has been tested by some early studies based on 
multijet + missing Et events originated by the associated production of the Higgs with a 
gauge boson or jet. This kind of studies put an upper bound on the invisible branching ratio 
of the Higgs of the order of 0.5 [1]. More recent studies are aiming to infer the viability of a 
non zero Br{h — > xx) by fitting the current data relying on Higgs searches. A study of this 
type can be found e.g. in Ref. [26] in which a fit of the current data has been performed 
by adding to the decay channels into visible states a branching ratio into two DM particles 
made free to vary. The outcome of such analysis is that BR(h — > xx) > 0.4 is excluded at 
the 95% confidence level with the best fit pointing towards values close to zero. 

The second set of general bounds applied in our analysis refers to the searches of 
FCNC processes, focussing on the three most important observables, namely the anomalous 
magnetic moment and the decay processes 6—7-57 an( i B s —> fJ, H~ ■ Regarding the 
anomalous moment of the muon, the SUSY contributions come from neutralino-slepton 
loops and chargino-sneutrino loops [27] with the latter typically dominating [28]. In absence 
of extra sources of flavor violation the most important parameter turns out to be the mass 
scale of the sleptons. Being this not directly related to DD we will not consider in the scans 
its experimental value but rather only an upper bound, together with the condition fi > 
which enforces the preferred sign of the SUSY contribution. 

Contributions to the b — > sj processes are instead induced by charged Higgs-top 
loops and chargino-stop loops [29] which have to be combined with the SM contribution, 
BrsM{b S7) — 3.0 x 10 -4 . While the charged Higgs-top loops always add constructively 
to the SM contribution, the behavior of chargino/stop loops depends on the sign of the 
quantity Re(fiAt). Having fixed the parameter to a positive value, they add destructively 
to the SM one when At < and constructively when At > . We also remark that the 
chargino-stop loops contribution is enhanced at high tan f3. 

Finally, regarding B s — > the additional contributions with respect to the SM 

arise from neutral heavy Higgs loops and are proportional to tan 3 /3 and to Atfi. This 
dependence on tan 3 j3 and the heavy Higgs mass makes, as already stated, this process 
also relevant for the Higgs phenomenology. In our analysis we made use of the numerical 
package Superlso [30, 31] for computing the rates of these three process imposing the range 
[2.77,4.37] x 10~ 4 for the values of Br(b -> 57) and upper bounds of 4.79 x 10~ 9 for 5a^ 
and 4.5 x 10~ 9 for Br(B s H + yT) [32, 33]. 

From the flat scans of the pMSSM parameter space performed varying the 12 para- 
meters previously introduced and accounting for all the mentioned experimental bounds, 
we selected four benchmark points for the Monte Carlo analysis on which our strategy to 
combine the DD technique and DM searches at the LHC is based. Before introducing these 
benchmark points, in the next section we will describe which DM configurations provide 
cosmologically viable scenarios in light of the recent DD limits from the XENON100 ex- 
periment (and current limits on SUSY searches) and how the DM spin-independent (SI) 
scattering cross section depends from the underlying parameters. The main quantities for 
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this study have been numerically computed through the package DarkSUSY [34]. 
2.2 DM relic density and spin-independent cross section 

The limit on the neutralino relic density is satisfied in rather different regimes, depending 
from whether the lightest neutralino is higgsino- and wino-like, or bino-like. In the former 
case, since higgsino and winos have unsuppressed couplings to W bosons via their charged 
counterparts, the pair annihilation rate is very large and the relic density matches the DM 
density in the Universe only going to the heavy mass regime, at about 1.4 and 2.7 TeV 
for, respectively, a pure higgsino and a pure wino [35]. Such regime (as much as for the 
couplings) is however not favorable for a DD analysis and will be discarded here. 

Moving to bino-like DM, we have an opposite situation, since bino-like DM couples to 
fermions only but the s-wave pair annihilation rate is helicity suppressed, namely propor- 
tional to the square of the mass of the final state fermion. As a consequence the value of 
the DM relic density exceeds the experimental limit with the exception of some peculiar 
regions of the parameter space. Interestingly these regions typically enforce definite rela- 
tions among the SUSY parameters which can be used in correlating an hypothetical DD 
signal with the collider detection. Some of these configurations are already well known in 
the context of CMSSM/SUGRA, thus for simplicity we will refer to them in the jargon 
borrowed from this class of models. 

Bino-like DM can be cosmologically viable in the following cases: 

• in presence of at least one rather low mass sfermion, namely with mass of the order of 
100-200 GeV, in such a way the annihilation cross section into fermion pairs matches 
the required value for a thermal WIMP. This is customarily dubbed bulk region. For 
the SUGRA/CMSSM it is ruled out by current LHC limits while in the pMSSM is 
still partially viable the possibility of light sleptons as opposed to higher mass squarks 
in agreement with the LHC bounds. 

• in the so called Higgs funnel region, in which there is an enhancement of the anni- 
hilation process, due to a s-channel resonance. In the CMSSM this can happen with 
the CP-odd Higgs, when m\ ~ 2m x . In the pMSSM it can also occur, at light DM 
masses, because of resonances associated to the Z boson and to the light CP-even 
Higgs. 

• in case the DM is a mixture of a bino and a small higgsino or wino component. The 
first case also occurs in the CMSSM in the so called focus point region while the 
pMSSM admits both the two configurations which are more generically dubbed as 
well-tempered neutralino [36]. 

• in the so called coannihilation scenarios, namely when one or more particles almost 
degenerate in mass with the DM contribute to keeping the DM in chemical equilibrium 
in the Early Universe. We also notice that coannihilating particles evade most of the 
current SUSY searches because the low mass splittings disfavor the detection of their 
decay products, after an eventual production at collider. 
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Figure 1: Curves at the same values of the DM bino (higgsino) fraction G g (Z g ), as re- 
ported in the plot, in the plane (M x ,asi)- This curves have been obtained by considering 
a simplified model featuring only M\ and \i as free parameters while all the remaining 
mass spectrum is fixed at about few TeV. In addition and tan/3 have been fixed to, 
respectively, 126 GeV and 10. The solid and the dashed black line are, respectively, the 
XENON100 limit and the XENON1T projected sensitivity. 



Moving to direct detection, looking at the general form for the neutralino-quark effective 
coupling in the small momentum transfer limit, reported e.g in [37] 4 , we can identify two 
main kinds of contributions determining the SI cross section. A first contribution depending 
on the SUSY parameters involved in the Higgs sector comes from the t-channel interaction 
mediated by a CP-even Higgs. In the decoupling regime only the light Higgs state h is 
relevant with the coupling scaling with the gaugino-higgsino mixing being zero for a pure 
bino or a pure higgsino. This effect can be seen in Fig. 1, where, in a simplified model 
featuring only Mi and fx as free parameters, we show the isocontour for the bino (higgsino) 
fraction G g (Z g ) in the plane DM mass versus SI cross section. 

The second kind of contribution comes from interactions mediated by squarks which is 
thus suppressed, compared to the former, by their mass scale. There is however a possible 
exception occurring in case of a very low mass splitting between a squark and the DM 
particle. Remarkably this setup suggests an interesting correlation with coannihilation 
scenarios. Moreover a suitable enhancement of the squark mediated interactions requires a 
sizable mixing between left and right-handed states. In our framework this condition can 
be fulfilled by third family squarks 5 . This is one relevant example of how the interplay 
of different observables may enforce correlations among SUSY parameters. On the other 

4 The expression there reported may not be strictly valid in case of coannihilations and need to be 
eventually corrected according to [38, 39]. 

5 This further motivates our choice of limiting to third generation squarks the possibility of being at low 
mass 
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Figure 2: The four benchmark points studied in this paper are shown together with the 
present Xenon exclusion limit and an estimation for the expected sensitivity of the Xenon 
experiment in its 1-ton configuration. 

hand the parameters of the stop mass matrix (related in our setup also to sbottom masses) 
determine as well the quantities Xt and M s entering in eq. (2.2). As already mentioned the 
constraint on the Higgs mass favors stops (and sbottoms as well) far from the coannihilation 
regime with the exception of a rather narrow range of values of \Xt/M s \ [21]. At the same 
time also the bounds form 6—7-57 seriously impact this scenario. 

2.3 Benchmark points for Direct Detection 

We now introduce four benchmark points selected for being representative of different cor- 
relation patterns between DD and collider observables. The first benchmark, labeled by 
xenl has been designed for the study of light DM scenarios and features a DM mass close 
to 56 GeV thus lying around the light CP-even Higgs resonance. We have then defined 
the point xen2, featuring DM-slepton coannihilations, with the aim of investigating the 
collider signals possibly induced by light sleptons. Finally, we have considered two bench- 
mark points, namely xen3 and xen4, corresponding to scenarios where DM is thermally 
produced through coannihilations with a light squark. The first of them, namely xen3, has 
been designed to have a very low value of the SI cross section and a mass close to the current 
exclusions limits, including the bounds from direct SUSY searches at LHC. The benchmark 
point xen4 is instead the heaviest benchmark considered in this paper, having a mass close 
to 500 GeV. The main features of the four benchmarks are summarized in Tab. 1 while 
Fig. 2 shows the four benchmarks in the plane DM mass versus SI cross section together 
with the present Xenon exclusion limit and an estimation for the expected sensitivity of 
the Xenon experiment in its 1-ton configuration. 
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m x 


m A 


A» 


m t\ 


"V, 


m- g 


mf 


&S I 


N 


xenl 


56 


890 


291 


557 


735 


815 


165 


2.3 x 10- 45 cm 2 


170 


xen2 


155 


1550 


353 


1715 


1777 


2254 


163 


3.12 x 10- 45 cm 2 


129 


xen3 


194 


2900 


1500 


340 


209 


2278 


732 


6.1 x 10- 46 cm 2 


24 


xen4 


477 


2780 


718 


527 


673 


665 


656 


1.35 x 10" 44 cm 2 


187 



Table 1: Summary table of the benchmark points. In the last column we report the 
number of recoil energies simulated in the DD analysis of section 3. 

Starting from these benchmark points we will illustrate how following our approach an 
hypothetical signal in the DD technique can be translated into a prediction for certain LHC 
observables. 

3 Direct detection 

In this section we will investigate the capability of a ton-scale direct detection experiment in 
reconstructing mass and SI scattering cross section associated with the pMSSM benchmarks 
previously introduced. Direct detection experiments seek to measure the nuclear recoil 
energy deposited in an underground detector when an incident % particle of the Milky Way 
DM halo interacts with a nucleus of the target material. In the case of spin-independent 
X-nucleus interactions the differential cross section for this process can be written as 

where Q is the energy transferred to the nucleus of mass Mjy from a \ particle moving 
with velocity u in the detector rest frame and A is the nucleus mass number. The form 
factor J~n(Q) accounts for the internal structure of the nucleus while a si is the x-nucleon 
SI scattering cross section at zero momentum transfer and M n = M x m p /(M X + m p ) the 
X-nucleon reduced mass. The differential rate of scattering events per unit detector mass in 
a direct detection experiment is given by the convolution of the incident flux of x particles 
with the differential cross section (3.1): 

dR _ a S i , 2 -2 ((D) f ,3,. Px-fxi^t) ( o 2) 

dQ-2M x Ml A:FN{Q) j H>v J U |u| ' (3 - 2) 

where the lower limit in the integral is the minimum velocity required for a \ particle 
to deposit the energy Q. In the case of an elastic scattering, such as for WIMP-nucleus 
interactions, the minimum velocity is given by 



QM N 

Vmin ~ V ~2Mf ■ {SJ) 

In Eq. (3.3) M r = M X M^ / [M x + Mn) is the x- nuc l eus reduced mass while M x represents 
the mass of the x particle. 
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The product p x ■ f x = of the local DM density p x and velocity distribution f x is 
equal to the phase-space density of the halo DM particles evaluated at the detector 
and in its reference frame. The local x phase-space density is computed assuming that 
the x particles account for the entire DM component and applying Eddington's inversion 
formula [40] to the mass model for the Milky Way introduced in section 3.1. This model is 
characterized by a spherically symmetric and isotropic distribution of DM particles and it 
furthermore assumes spherical symmetry for the underlying Galactic gravitational potential 
<£(r). In this limit, Eddington's formula [40] gives an one-to-one correspondence between 
the DM halo density profile Ph{ r ) and its phase-space density F^: 



£ d 2 p h (W 1 fd Ph 



(3.4) 



where we introduced the relative potential 'l'(r) = — $(r) + $(r = R v i r ) and the virial radius 
Rvir, corresponding to the radius at which the DM halo is truncated. Under the assumptions 
considered here, F^ is a function of the relative energy £ = —E+Q(r = R v i r ) = — E\ i \ n +'fy(r) 
only, with E and E^ n , respectively, the total and kinetic energies. The relative energy 
depends on the DM particle velocity and position in the Galaxy and therefore determines the 
space-time dependence of Fh itself. To numerically compute Fh at the local Galactocentric 
distance it is actually simpler to implement Eq. (3.4) by changing the integration variable 
from ^ to the radius of the spherical system r. Then, the radial dependence of ph and $ 
from the local position in the Galaxy out to R v i r is the only information needed to calculate 
Fh- Further details on the implementation and validity of this approach can be found in [4]. 
In summary, Eddington's inversion formula allows us to express the integrand of Eq. (3.2) 
directly as 

p x -f x (u,t) = F h [£ (v)} , (3.5) 

with v(t) = v bs(i) +u, being v Q b s the velocity of the detector in the Galactic frame, where 
Fh is computed, and £q the relative energy at the local Galactocentric distance. For v G b s 
we take the expression 

Vobs(t) = v LS r + v 0>pcc + V©(t) , (3.6) 

where the Local Standard of Rest velocity vlsr (namely the rotation of the Sun with 
velocity 0o around the Galactic center; calculable given a mass model for the Galaxy), the 
Sun's peculiar velocity VQ )Pec and the motion of the Earth around the Sun Vq are given 
by 

v LS R = (0,G ,0), 

v 0,pec = (U Q ,V Q ,W Q )~ (10,5.2,7.2) km/s 

V e (t) = V( B [iicosL}(t-t 1 ) + e 2 smuj(t-t 1 )] , (3.7) 

with V® = 29.8 km/s and t\ = 0.218 is the fraction of the year before the Spring equinox. 
Finally, i\ and £2 are the directions of the Earth's velocity at times t\ and t\ + 0.25 years. 
In Eq. (3.7) all velocities are expressed in Galactic coordinates, where x is the direction 
to the Galactic Center, y the direction of disk rotation and z identifies the North Galactic 
Pole. 
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3.1 Galactic model 

The Galactic mass model adopted in the present analysis to compute Fh has been ex- 
tensively investigated in Refs. [3, 4]. For completeness, we briefly summarize it in what 
follows. The model consists in two luminous mass components, namely the stellar disk and 
the Galactic bulge, and in a DM halo. The Galactic model parameters introduced in this 
section will be then treated as nuisance parameters in the Bayesian analysis of our pMSSM 
benchmarks performed in the next subsections. 

Regarding the stellar disk, we assume a mass density profile which in cylindrical coor- 
dinates (R, z) with origin in the Galactic center is given by 



p d (R, z) = —^e R d sech 2 ( — 
2z d \z d 



with R < R, 



dm 



(3.8) 



where is the central disk surface density, R d and z d are length scales in the radial 
and vertical directions, while R dm is the truncation radius of the disk. R d m is assumed 
to scale with the local Galactocentric distance Rq according to the prescription R dm = 
12[1 + 0.07(i?o — 8 kpc)] kpc and the vertical scale z d is fixed to the best fit value suggested 
in Rcf. [41], z d = 0.340 kpc. 

The bulge/bar region is instead characterized by the mass density profile: 



p bb (x,y,z) = p bb 



-1.85 



exp(-s a ) +exp 



where 



q 2 b (x 2 + y 2 ) + z 2 



and 



+ 



+ 



(3.9) 



(3.10) 



We implement in this analysis an axisymmetrized version of Eq. (3.9), and assume Xb — 
yb = 0.9 kpc- (8 kpc/i?o), z b = 0.4 kpc- (8 kpc/i?o) and qb = 0.6. See also [3] concerning the 
choice of these parameters. Rather than using the two mass normalization scales T, d and 
Pbb as free parameters, we re-parameterize these in terms of two dimensionless quantities, 
namely, the fraction of collapsed baryons /b and the ratio between the bulge/bar and disk 
masses V: 

n BM + n b M bb + M d + m Hi + M Ha 



fh 



Mbb 
M d 



(3.11) 
(3.12) 



In Eq. (3.11) we also included the sub- leading contributions to the total virial mass M v i r 
(defined in the following) associated with the atomic (Hi) and the molecular (H2) Galactic 
gas layers, with profiles as given in [42]. In summary, the free parameters describing the 
luminous components are Rq, R d , fb and T. 

Concerning the DM halo component we consider an Einasto profile [43, 44], which is 
favored by the latest N-body simulations and is given by 



Ph(r) = p'fE {r/a h ) , 



(3.13) 
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with 

f E (x) = exp 



— {x aE - 1) 

OLE 



(3.14) 



The reference normalization p' and the scale radius ah in Eq. (3.13) are often rewritten as 
a function of the virial mass M v { r and of the concentration parameter c v i r by inverting the 
relations: 

i «■ 47T _ 3 OpM + ^ b , /' Rmr t 2 / \ /q -i r\ 

-Mmr = -5- A„ ir p -K^r- = o 4?r / « rr /M r ) ( 3 - 15 ) 

3 "DM Jo 

c„ ir . = R vir /r^ 2 (3.16) 

where the virial overdensity A„j r in the first equation is computed according to Ref. [45] 
while po is the mean background density today. The presence in this equation of £1dm 
and f^b, the DM and baryon energy densities in units of the critical density, reflects our 
assumption that only a fraction equal to ^dm/(^dm + ^b) of the total virial mass consists 
of DM. Their values have been set according to the mean values from the fit of the 7- year 
WMAP data [46]. In the second equation, instead, r_2 is the radius at which the effective 
logarithmic slope of the DM profile is equal to —2. Finally, we assume that the baryons 
which do not collapse in the disk are distributed according to the same profile of the DM 
component. The free parameters describing the DM halo are therefore M v i r , c v i r and a. 

The Galactic model used to compute Fh also includes an additional parameter, the 
so-called anisotropy parameter /3 = 1 — of /of , where a r and at are respectively the radial 
and tangential velocity dispersions of the halo stars. It has been introduced in Ref. [3] to 
include in the parameter estimation the velocity dispersion measurements of Ref. [47] . 

3.2 XENON1T 

In a real experiment Eq. (3.2) has to be modified in order to account for experimental 
limitations related to the efficiency of the detector, its finite energy resolution and the 
presence of an experimental noise. In our discussion of the capability of a ton-scale direct 
detection experiment in reconstructing the DM mass and SI scattering cross section we will 
implement the experimental limitations which would apply, as a first approximation, to the 
XEN0N1T experiment. We will therefore consider a constant efficiency e = 0.3 and include 
finite energy resolution effects by assuming a Gaussian probability 

Q) = 1 e -{E-Q?/2„UQ) (3 17) 

V 27r 4e(Q) 

that to an event of energy Q is wrongly assigned an energy in the interval E + dE. In 
Eq. (3.17) the energy dependent dispersion <7xe(Q) is the one estimated for the Xenon 
detector [48]: 

a X e(Q) = y / Q/KeV(0.579 KeV) + 0.021Q . (3.18) 
The observed direct detection rate is hence equal to 

jn poo jo 
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Regarding the experimental noise, we consider two contributions characterized by a different 
energy dependence: a flat noise component (i.e. constant in energy) and an exponential 
noise component oc exp(— E/Eq) with Eq = 10 KeV. Finally, in all the conclusions drawn 
below we will assume an experimental raw exposure of 1000 kg-years and a signal region 
ranging between 8 KeV and 75 KeV. 

3.3 Bayesian forecasts 

To test the capability in reconstructing the DM mass and SI cross section of a ton-scale di- 
rect detection experiment with specifications resembling the ones expected for XENON1T, 
we will first simulate a set of mock direct detection data and then analyze them employ- 
ing a Bayesian approach supported by a Markov Chain Monte Carlo (MCMC) Likelihood 
scanning technique. Mock data have been simulated for each pMSSM benchmark analyzed 
in this work. 

As already mentioned, in the context of the direct detection we need to extend the 
pMSSM parameter space in order to include additional parameters crucial for the compu- 
tation of the differential rate (3.2). These are the parameters introduced in section 3.1 
when defying the Galactic model used to compute the DM phase space density Fh within 
Eddington's approach. We denote by p the combination of the pMSSM parameters used 
to describe a certain benchmark and the parameters characterizing the underlying Galac- 
tic model. Similarly, by p^, we denote the value of p at the considered benchmark point. 
Whereas the assumed pMSSM parameters differ for the four benchmarks studied in this 
paper, in all cases we adopt the same benchmark Galactic model. This has parameters 
fixed to the mean values found in [4]. The mean local DM phase space density for the 
Einasto, NFW and Burkert profiles can be downloaded from the online material associated 
with Ref. [4]. 

The data simulation consists in sampling a set of N recoil energies Ei, i = 1, . . . , N 
from the probability density function f(E, p&) obtained evaluating at the benchmark point 
Pb the spectral function 

HE r>)= 1 ( dR(E,p) 1 exp(-E/E ) 

Mp) I dE b-a E [exp(-a/E Q ) - exp(-6/£ )] 

where a =8 KeV, b =75 KeV and //(p) = J^dEdR(E,p)/dE + 2. This choice corresponds 
to assume two background events in the signal region during the exposure time: 6 one from 
the flat noise component and one from the exponential noise component. For a given point 
p in parameter space, Eq. (3.20) represents the probability that the XENON1T experiment 
measures an energy (infinitesimally close to) E within the signal region (a, b). The first term 
on the right hand side of this expression corresponds to the contribution to this probability 
from the true DM signal while the second and the third terms are associated with the 
flat and exponential noises respectively. The number of recoil energies N is drawn from a 

6 As expected for ton-scale direct detection experiments which aim at reducing the number of background 
events measured during the exposure time to ~ 1. 



(3.20) 
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Poisson distribution of mean A& with probability 

Poisson(iV; A 6 ) = exp (-A 6 ) (3.21) 

where = /^(pb) is the number of events (i.e measured recoil energies) expected for the 
benchmark under analysis. In Tab. 1 we report the values of N generated in the simulations 
analyzed in this work. 

Having explained how to generate mock data from a given benchmark point, we can 
now focus on their Bayesian analysis. As already mentioned the aim is to reconstruct the 
dark mater mass and SI cross section. We will estimate the "preferred values" for these 
quantities by determining their marginal posterior probability density functions (PDF). 
According to Bayes' theorem the posterior PDF of certain parameters p is proportional to 
the product of the Likelihood function £(p, d) and the prior probability density 7r(p): 

P(p , d) ^£^Pl 0.22) 

where d is the array of datasets used to constrain the parameter space. The Bayesian 
evidence £ (d), being independent from p, plays the role of a normalization constant when 
performing parameter inference. The marginal poster PDF of a generic function g of the 
parameters p, e.g. a$i and M x in our case, is given by the expression 

p(g\d) = J dp %(p) - g) P(p|d) , (3.23) 

which follows from the definition of conditional probability. MCMC scanning techniques 
have been then used to sample from Eqs. (3.22) and (3.23) the marginal posterior PDF of 
the DM mass a SI cross section. 

The Likelihood function is a key ingredient of the present analysis and it has been 
constructed following the prescriptions given in Ref . [48] . The contribution to the Likelihood 
from the simulated direct detection data is given by 

- log£ dd (^,iV,p) = -N- iVlog (^M) - f>g (y||^y) (3-24) 

and it has been normalized in such a way that logCdd(Ei, N,p) = 0, where /x(p) = N. 
This expression reflects both the information on the number of simulated recoil events N, 
included in the Poissonian pre-factor, and the information about their energies Ei, encoded 
in the spectral function (3.20). In addition to the direct detection contribution, we have 
included in the Likelihood a term which assigns a higher statistical weight to points in 
parameter space characterized by a "physical mass spectrum". This term implements in 
the analysis the LEP2 constraints on SUSY parameters, bounds from collider searches 
of Higgs decays into r pairs (which can be interpreted as an mA-depedent upper bound 
on the value of tan/3) and, finally, a lower bound of 600 GeV on the gluino mass [23]. 
These constraints, being expressed in terms of lower limits £ r (d) on certain parameters 
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or theoretical predictions r r (p), generate a contribution to the Likelihood function of the 
following form 



-lGg£«(d,p) = 



where the index r labels the constraints on the pMSSM mass spectrum previously mentioned 
and o>(d) is an estimate of the error associated with the value of the lower limit £ r (d). 
We also included in the Likelihood a term accounting for the recently discovered boson 
compatible with the Standard Model Higgs: 

-log/:m gg s(d,p) = - j (3.26) 

where m/ l (p) is the pMSSM prediction for the light CP-even Higgs boson, rhh(d) = 126 
GeV and tr^(d) = 3 GeV. Summing up all the contributions, the final form of the Likelihood 
function considered in this analysis is 

- log £(q, p) = - log C dd (Ei ,N,p)-J2 log 4c ( d > P) " lo § AWd, p) . (3.27) 

r 

The last information needed to complete the analysis setup concerns with the choice 
of the prior PDF 7r(p). In the literature the most investigated choices correspond to flat 
and logarithmic priors PDF in parameter space. Regarding the pMSSM mass parameters, 
in the reconstruction of the DM mass and SI cross section we will consider both cases 
separately to assess the impact of priors on the final results. Concerning instead tan f3 
and the trilinear couplings, following [49], only flat priors will be considered. To carefully 
sample the posterior PDF of the pMSSM parameters, we extended in the present analysis 
the ranges introduced in Eq. (2.1), varying squark and slepton masses between 50 GeV and 
4 TeV. The ranges considered for the other parameters coincide instead with the ones of 
Eq. (2.1). Galactic model parameters will be treated as nuisance parameters, using the 
knowledge on their distributions acquired in Refs. [3, 4]: the posterior PDF found in these 
references for the Galactic model parameters through a compilation of various dynamical 
constraints will be considered in the present analysis as the prior PDF of the astrophysical 
parameters. As a final remark, we mention here that all the parameters controlling the 
hadronic matrix elements required to evaluate the DM-nucleus SI scattering cross section 
have been fixed at the mean values of Ref. [48]. In fact, given the high dimensionality of 
the pMSSM and the limited constraining power of the data, we preferred to consider in this 
paper as nuisance parameters only the astrophysical parameters to which we could assign 
prior PDF motived by previous works. We refer the reader to Refs. [50] for an analysis 
assessing the role of hadronic matrix elements in the context of DM direct detection. 



1 f Tr(p)-e r (d) ^ ' 

oy(d) 



if r r (p) < £ r (d) 

(3.25) 

otherwise 
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Figure 3: Histograms showing the distributions of the best fit DM masses derived by 
fitting 1000 data realizations with the theoretical expectation (3.20). The left panel refers 
to benchmark xenl and the right panel to benchmark xen2. The red cross represents 
the value of the DM mass at the benchmark point. 



3.4 Statistical limitations 

By reconstructing the DM mass and SI cross section we mean what follows: to sample 
the posterior PDF (3.22) through a MCMC algorithm 7 and derive from it the marginal 
PDF of the DM mass and SI scattering cross section. We will present results in terms of 
a% Bayesian credible regions containing therefore a fraction a of the posterior probability. 
Focusing on credible regions rather than on frequentist confidence intervals allows to cir- 
cumvent possible issues regarding the coverage of the resulting intervals, a problem which 
has recently attracted a great attention in the context of the Minimal Supersymmetric 
Standard Model [51] and in the DM direct detection field [11]. 

There is however another possible source of "statistical limitations" affecting the re- 
construction procedure described above. This concerns with the fact that the mock data 
implemented in the analysis, for instance the recoil energies Ei, represent just one single 
realization of the true underlying distribution of the data. If the simulated data do not pro- 
vide a faithful representation of their true distribution, because of statistical fluctuations, 
the resulting reconstruction can be significantly biased. In Ref. [11] this effect has been 
observed for DM candidates covering a non negligible spectrum of masses and SI scattering 
cross sections. It is therefore important to check the sensitivity of our benchmark recon- 
structions to the particular realization of the data implemented in the analysis. Given the 
sophistication of the Galactic model of section 3.1 and of the particle physics parameter 
space considered in this work, to run several MCMC for different realizations of the data is 
not an efficient way of achieving this task. Therefore, to asses the dependence of our results 
from the use of a single data sample for each benchmark we followed a different approach 
based on a simple one dimensional x 2 minimization: we first sampled for each benchmark 
a large number of datasets (N,Ei) as explained in the previous subsection. Then, for each 

7 In this analysis we employ a modified version of the SuperBayes code interfaced with DarkSUSYroutines 
to compute direct detection observables. 
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Figure 4: As in Fig. 3 but the left panel refers to benchmark xen3 and the right panel to 
benchmark xen4. 



realization of the direct detection data, we fit the resulting distribution of the ./V energies 
Ei against the theoretical expectation (3.20) varying the DM mass only and fixing the SI 
scattering cross section and astrophysical parameters to their benchmark values. The re- 
sulting distribution of best fit values for the DM mass will allow us to estimate the degree 
of sensitivity of the benchmarks studied in this paper to the specific data sample considered 
in the reconstruction of the DM mass and SI cross section. It is in fact the DM mass the 
quantity mostly affected by statistical fluctuations in the distribution of the simulated recoil 
energies. 

The results of this analysis are shown in Fig. 3 for benchmarks xenl and xen2, while 
Fig. 4 refers to benchmarks xen3 and xen4. These histograms show the distributions of 
the best fit DM masses derived through the fitting procedure explained above for 1000 
different data realizations. In each figure the red cross represents the value of the mass at 
the benchmark point. In all cases we see that the distribution is not perfectly symmetric 
around the benchmark point: there is a general trend in overestimating the DM mass. 
This result is in agreement with what found in Ref. [11], where the impact of statistical 
limitations on the reconstruction of the DM mass and SI cross section has been studied 
for a broad spectrum of DM candidates. This study shows that the most sensitive DM 
candidates are those with high masses and low SI cross section, i.e. candidates associated 
with a relatively low number of expected direct detection signal events. This trend is 
confirmed by the outcome of our analysis, which shows how the benchmarks potentially 
more sensitive to statistical limitations are those with the highest mass, namely xen3 and 
xen4, as clearly reported in Fig. 4. The reason behind this phenomenon is that when the 
number of recoil events is not sufficiently large to produce a faithful description of the 
underlying recoil energy distribution, a small excess of events at high energy, resulting from 
a statistical fluctuation, is sufficient to alter the result of the fit: this excess of events at 
high energy will be better fitted by an artificially large value of DM mass (see also Fig.l in 
Ref. [11]). 

For the benchmarks studied in this paper we notice however that this effect, though 
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Figure 5: Two-dimensional marginal posterior PDF in the plane DM mass versus SI 
scattering cross section. The left panel refers to xenl while the right one to xen2. In 
each panel the colored bar follows the posterior PDF, arbitrarily normalized to one, from 
its peak down to zero. Black contours correspond to the 68% and 95% credible regions 
obtained implementing in the analysis flat priors for the pMSSM mass parameters, while 
the green contours correspond to the analogous quantities derived from log-priors. The 
white dot and the blue cross represent respectively the posterior mean (from log-priors) 
and the benchmark point. 



present, it does not influence dramatically the best fit distributions of the DM mass. Since 
therefore we do not observe any strong pathology related to the selected benchmarks, 
we continue our analysis focusing on their mass and SI cross section reconstruction at 
XENON1T. 

3.5 Mass and cross-section reconstruction 

The procedure to reconstruct the DM mass and SI cross section introduced in the previous 
subsections has been applied to the benchmarks xenl, xen2, xen3 and xen4. For each 
benchmark we focused on a single data realization considering both flat and log-priors 
for the pMSSM mass parameters. In Fig. 5 and 6 we show for these benchmarks the 
two-dimensional marginal posterior PDF in the plane DM mass versus SI scattering cross 
section. In each panel the colored bar follows the two-dimensional posterior PDF which 
has been conveniently normalized to one. Black contours correspond to the 68% and 95% 
credible regions obtained implementing in the analysis flat priors for the pMSSM mass 
parameters, while the green contours correspond to the analogous quantities derived from 
log-priors. Finally, the white dot and the blue cross represent respectively the posterior 
mean (from log-priors) and the benchmark point. 

The plots confirm the general tendency in overestimating the DM mass described in 
the previous subsection, here evidenced by the high mass tail featured by the contours and 
partially by the mismatch between the posterior means and the values of the mass and 
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Figure 6: As in Fig. 5 but the left panel refers to benchmark xen3 and the right panel to 
benchmark xen4. 



SI cross section of the original benchmark. On the other hand we observe a significant 
improvement of the reconstruction performance when adopting log-priors in the MCMC 
scanning procedure. Indeed, this choice tends to "compensate" possible limitations in the 
DM mass reconstruction, due to statistical fluctuations in the distribution of the simulated 
data, by giving more statistical weight to low values of the parameters [52] (and therefore 
to the DM mass), thus partially reducing the high mass tail of the marginal posterior PDF. 
Because of this compensation mechanism we argue that log-priors provide in this context a 
more faithful description of the actual potential of upcoming 1-ton DD experiments and we 
will therefore employ results obtained with this choice when translating our DD analysis 
into a forecast for the LHC. 

For a given benchmark, in addition to this effect, the quality of the reconstruction is 
also affected by the value of the DM mass itself. Indeed, the differential recoil rate, for DM 
masses largely exceeding the one of the nuclei in the target material, depends on the DM 
parameters only through their ratio osi/M x . Clearly, this degeneracy produces credible 
regions spanning over "erroneously" large portions of the plane (M x ,asi)- In the present 
analysis the benchmarks most largely affected by this degeneracy are benchmarks xen3 
and xen4, as also shown in Fig. 6. 

In the next section we will translate the preferred regions in the plane DM mass versus 
SI cross section determined here into a prediction for certain classes of signals at the LHC. 
Since most of the relevant SUSY production cross sections are very steep functions of the 
pMSSM mass spectrum, and therefore of the DM mass, also the LHC analysis described 
in the following will be indirectly affected by the statistical limitations producing large 
uncertainties in the DM mass reconstruction. Large exposures are therefore crucial to 
enhance the effectiveness of the strategy illustrated here to combine the direct detection 
technique with DM searches at the LHC. 
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In summary, a ton-scale experiment resembling the XEN0N1T specifications will be 
able to achieve a robust determination of mass and SI cross section for DM candidates 
similar to the one associated with the benchmark point xenl. In fact, in light of our 
simulations, xenl clearly emerges as the most promising benchmark considered in this 
work. This result follows from the large number of recoil events in the simulation of xenl 
(170) and from the fact that light candidates are in this context simpler to identify with the 
direct detection technique, having a steeper energy spectrum which minimizes the previously 
mentioned statistical limitations. The benchmark points xen2, xen3 and xen4 feature 
instead a less efficient reconstruction but still in moderate agreement with the true values 
of the DM parameters related to the original benchmark. These points are in fact more 
sensitive to the already mentioned statistical limitations and to the mass/cross section 
degeneracy, typical of heavy DM candidates. 

4 Collider analysis 

Having derived for each benchmark the accuracy within which the relevant DM properties 
are reconstructed by a 1-ton DD experiment resembling XENON1T, we now focus on how 
to translate this information into a prediction for certain LHC observables. Our analysis 
consists in two steps: first we will reprocess the Markov chains produced in the DD analysis 
applying to the models in the chains additional constraints from the relic density, limits 
on the muon anomalous magnetic moment as well as bounds from B-physics processes. 
Then, when possible, we will determine the production processes most likely associated to 
a subsample of points in the chains passing these additional constraints and study the cor- 
responding LHC response in its 7 TeV center of mass energy configuration. This procedure 
is the essence of the method proposed in this paper to correlate the DD technique with DM 
searches at the LHC. 

To rigorously assess the impact on the DD analysis of section 3 of additional observables 
one should run new MCMC scans with an extended Likelihood properly incorporating the 
new constraints. Although this is by now technically achievable in the most simple SUSY 
realizations (see [52, 53] for some examples) this procedure would be computationally trou- 
blesome to implement in the present analysis because of the need of repeated lengthy 
computations of quantities like the DM relic density, especially in the context of a high di- 
mensional MSSM parameter space, as the one considered in this work, and when studying 
points in parameter space achieving the correct relict density via co-annihilation processes. 
For these reasons we have employed the less computationally demanding procedure of ap- 
plying a posteriori the limits from the relic density and flavor processes to the Markov 
chains generated in the analysis of section 3. More precisely, for each benchmark we have 
computed the observables just mentioned only for the PDF samples enclosed in the 95% 
log- prior contours reported in Figs. (5) and (6). The limits implemented for <Ja„, b — > s"f 
and B s — > are the same described in section 2. Regarding the relic density, instead, 

we have just required that it does not exceed the cosmological limit. This is motivated by 
the fact that too strict requirements may excessively affect set of points not informative on 
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Table 2: Summary table of the SUSY searches employed in our collider analysis, 
this observable 8 . 

As already mentioned, the additional constraints imposed on the model parameters 
restrict the 95% credible regions determined in section 3 and produce the subregions on 
which our collider analysis has been mostly focused. We have then identified the production 
processes most likely associated to the SUSY models contained within these subregions and 
employed accordingly the analysis strategies applied by the ATLAS [59] collaboration to a 
few points sampled from them. Indeed, for each of these points we have simulated the signal 
possibly detectable by the ATLAS experiment and obtained the number of expected signal 
events for several combinations of final states, as summarized in tab. 2, at the luminosities 
relative to the current available data. From the comparison with the current limits we have 
determined whether the various configurations are already severely constrained or will be 
probed by the very next upgrades in luminosity. 

The signal events have been generated through the numerical package Madgraph [60] 
combined with the detector simulator PGS normalizing the production cross-section at 
the NLO through the package PROSPINO [61]. The simulations have been validated by 
generating the main backgrounds for the considered search channels and cross-checking our 
results with the number of background events inferred by the data. The number of events 
obtained by the simulations has resulted compatible with expectations. 

4.1 Analysis of xenl 

We start our analysis with the lightest DM mass scenario, namely xenl. Fig. 7, obtained 
selecting points in our Markov chains passing flavor and relic density constraints, evidences 
that values of the DM relic density compatible with current limits can be achieved in two 
setups only: the first is represented by the two narrow regions lying at values of the DM 
mass of around 45-60 GeV, thus coinciding with the original benchmark scenario. In this 
regions a bino-like DM candidate is driven to the correct relic density by resonance enhanced 

8 Alternatively one could also consider deviations from the standard WIMP paradigm, as enforced by 
non thermal production [54-56] or alternative cosmologies [57, 58]. This would require however a different 
treatment of the local DM density [53]. However, SUSY parameters are less directly correlated in these 
scenarios. 
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Figure 7: Points of the Markov chains associated to the benchmark xenl passing flavor 
and relic density constraints in the plane (M x ,Qh 2 ) following a color pattern related to 
the DM higgsino fraction. 



interactions mediated either by the light CP-even Higgs boson or, for slightly lower mass 
values, by the Z boson. In the second setup, instead, occurring at higher values of the DM 
mass and corresponding to the yellow region in the plot, DM is mainly a pure higgsino 
characterized by very low values of the relic density and by a value of the mass larger than 
100 GeV as required by the LEP bound. 9 

Notably, there is a clear correspondence between the second setup, represented by the 
yellow band in Fig. 7, and the high mass tail of the two-dimensional marginal posterior PDF 
in the plane DM mass versus SI scattering cross section shown in Fig. 5. As argued in section 
3.4, this tail arises as a consequence of statistical fluctuations in the recoil energy spectrum 
affecting the accuracy of the DM mass reconstruction. We therefore conclude that the 
higgsino DM setup can be effectively disentangled from the bino DM case (corresponding to 
the original benchmark) by improving the performance of the DD based parameter inference, 
for instance, by increasing the experimental exposure. In view of this consideration we have 
limited our collider analysis to a subset of models sampled from the PDF of section 3 and 
corresponding to gaugino DM. 

For light DM candidates, as the one considered in this section, the most relevant para- 
meters entering the SI cross section and relic density calculation are the DM mass itself and 
the value of the \x parameter. This implies that, from the collider point of view, the most 
interesting way to probe this scenario consists in the search for direct gaugino production, 

9 Fig. 7 also reports a very small amount of points at gaugino DM in the mass interval of around 70-300 
GeV. For these points a correct value for the relic density is guaranteed by either a low mass slepton or by 
a CP-odd Higgs. In the first case, however, the slepton should lie very close to the LEP limit and it is thus 
statistically disfavored by our Monte Carlo scanning procedure, while in the second low values of rriA are 
disfavored by Br(B 3 — s- 
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Table 3: Values of the non-SM visible cross section associated to 3 leptons + E™ 1SS 
search channel for 5 sample points relative to the light DM scenario depicted by the 
benchmark xenl. The results rely on two signal regions associated, respectively, to the 
decay of the gauginos into off-shell and on-shell gauge bosons. The experimental limits 
are a° b s s = 3fb and < s s = 2fb. 

namely the associated production of x% an d Xi • Assuming that there are no intermediate 
sleptons between these two states, they can decay either into Z, with the bosons which 
may be also off-shell, or in the pair Wh, in addition to missing energy. Our study has 
been focused on the first possibility employing the ATLAS dedicated search through events 
with 3 leptons in the final state [62]. A strategy for the detection of the Wh channel has 
been instead proposed in [71] but it requires much higher luminosities for being detected 
compared to the previous channel. 

Processes involving the production of colored particles do not play a relevant role, at 
least up to moderate values of M?, since the current collider limits imply heavy masses 
such that the relative production rates can be regarded as subdominant. Notice that this, 
in the pMSSM, is strictly valid for the gluino and the squarks of the first two generations; 
third family squarks can feature much less stringent limits 10 . This is likely not the case 
for such low values of the DM mass. Indeed sbottom masses are excluded up to around 
390 GeV at 2.051b" 1 for DM masses less than 60 GeV [75] and preliminary results from 
the 4.7fb _1 update [68] move the lower limit up to around 500 GeV. Regarding the stop 
(see, for reference, [65, 66, 76]), some values in the low mass range are still not excluded 
(the reason is that the stop features several different decay channels and then the study of 
very different signals, each one probing a limited range of masses, is required. More details 
will be provided in the following sections). However, for the values of the SI cross section 
emerging in the present scenario, squark mediated interactions are not expected to play an 
important role, then it is more interesting to focus on collider searches probing instead the 
relevant parameters, as for instance the previously mentioned 3 lepton final state. 

Tab. 3 reports the results relative to 5 models, named xenl-l,. . . , xenl-5, sampled 
from the PDF relative to xenl and passing relic density and flavor constraints. The table 
reports, for each model, the values of M2 and [i which determine the characteristics of 
the produced neutralino and chargino, together with the visible cross-sections for non-SM 

10 A part from the cases discussed in the next sections, in which it is established a correlation with the 
relic density requirement, we exclude as well the possibility of very compressed spectra (see e.g. [72-74] 
and references therein for some studies). 
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Figure 8: Missing energy distribution for the events (signal+background) including 
three leptons in the final state expected for the models xenl-1,. . . ,xenl-5 sampled from 
the chains relative to benchmark xenl. This result has been obtained in a run of the 
LHC with an energy in the center of mass of 7 TeV and a luminosity of 4.7 fb _1 . 

processes, i.e. the product ere A of the theoretical production cross-section a, the detector 
efficiency e and the acceptance A which accounts for the reduction of the signal strength 
due to the analysis cuts. DD combined with relic density and flavor requirements has a 
good constraining power for the [i parameter, which we have found to vary on a moderately 
restricted range with a lower bound of around 200 GeV, as expected from the functional 
form of the SI cross section, while leaves substantially M2 free to vary. The results in 
the table shows that the W Z channel can efficiently probe models with moderate values 
of M2 with some regions already excluded and some other which can be detected in the 
near future with an O(l) increase of the luminosity. At high values of M2, instead, the 
chargino and the second neutralino increase their higgsino fraction favoring the decay into 
Wh and requiring therefore a different search strategy compared to the one considered here 
involving three leptons in the final state. 

Fig. 8 shows the missing energy distribution corresponding to events (signal+ back- 
ground) including three leptons in the final state simulated from the models xenl-1,. . . , 
xenl-5 in a run of the LHC with an energy in the center of mass of 7 TeV and a luminos- 
ity of 4.7 fb _1 . In summary, the benchmark scenario xenl features a clearly identifiable 
signature in the three lepton final state emerging as a genuine forecast of our DD analysis 
which in this case can be effectively translated into the results of Fig. 8 and Tab. 3. 

4.2 Analysis of xen2 

The benchmark point xen2, similarly to xenl, has been designed to study MSSM realiza- 
tions featuring a light EW interacting sector (see e.g. [77] for a similar approach) evading 
current constraints which would instead apply to models including colored sectors lying in 
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Figure 9: Points of the Markov chains associated to the benchmark xen2 passing flavor 
and relic density constraints in the plane (M x ,Qh 2 ) following a color pattern related to 
the DM higgsino fraction. 



the multi-TeV region. As for the case of xenl, in a preliminary investigation of the bench- 
mark xen2, we selected from our Markov chains the models passing flavor and relic density 
constraints. These points are shown in Fig. 9 where again the yellow band corresponds to 
higgsino DM configurations. Contrary to the case of xenl, however, in the present scenario 
the higgsino DM models do not only populate the tail of the marginal posterior PDF in 
the plane DM mass versus SI scattering cross section: comparing in fact Fig. 5 with Fig. 9 
one can clearly observe that there are regions in the plane (M x ,asi) with a significant 
statistical weight corresponding to DM masses above 100 GeV and therefore representing 
a large fraction of points in the yellow band of Fig. 9. This is obviously related to the fact 
that the LEP bound implemented in the Likelihood of our MCMC scans is only effective 
for exotic states lighter the 100 GeV, while the DM candidate associated with xen2 has a 
mass of approximately 155 GeV. 

For these reasons the outcome of our DD based analysis cannot resolve alone the de- 
generacy among pure gaugino and pure higgsino DM configurations (see Fig. 9). A possible 
strategy to alleviate this problem could be to encode in the Likelihood the requirement of a 
DM relic density lying within 2-3 a from the cosmological limit; this requirement cannot be 
fulfilled by higgsinos with relatively low values of the mass, as required by this benchmark 
scenario. A potentially interesting alternative is, on the other hand, adding information 
coming from indirect detection. Indeed bino-like and higgsino DM feature different sig- 
natures from the point of view of indirect detection. In particular, current experimental 
facilities have reached enough sensitivity to probe light higgsino DM scenarios with already 
strong constraints coming form antiproton data [78] and new stronger bounds eventually 
arising in case that the claim of detection of a gamma-ray line in the Fermi-LAT data [79] 
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Table 4: Values of the non-SM visible cross sections associated to the 3 leptons + E™ 
search channel for 5 sample points relative to the DM scenario depicted by the bench- 
mark xen2. The results rely on three of the signal regions considered in Ref. [62]. The 
experimental limits are, respectively, u°^ s = 3fb and cr°^ s = 2fb. 
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Table 5: Values of the non-SM visible cross sections associated to the 2 leptons + E™ 1SS 
search channel for 5 sample points relative to the DM scenario depicted by the bench- 
mark xen2. The results rely on three of the signal regions considered in Ref. [63]. The 
experimental limits are, respectively, <7°^ s = 2fb, <t°^ s = 11.4 fb and cr°^ s = 2fb. 

will be confirmed (for more details see e.g. [80]). 

Leaving these possibilities for future works, we will focus in this section on the scenario 
corresponding to gaugino DM while a discussion on possible peculiar signatures of higgsino 
DM at the LHC is postponed to the next section. Moreover, we will restrict in the present 
collider analysis to the DM mass range 135-175 GeV, corresponding to a region of high 
posterior PDF in the plane DM mass versus scattering cross section able to capture the 
main features of the present scenario. Within this mass range there are several different 
possibilities to achieve a DM relic density compatible with the cosmological limit. The main 
ones are the coannihilation with a slepton and resonances occuring in annihilations mediated 
by the CP-odd Higgs, with the latter challenged by the constraints form B s — > /i + fi~ . In 
addition to these, we have also found configurations in which the parameters M\ and M2 
are very close to each other, featuring the already mentioned setup dubbed "well tempered 
neutralino": in this case the correct relic density is achieved by the interplay of a non 
negligible wino fraction of the DM and its coannihilations with the lightest chargino and 
the second neutralino. 

This variety of scenarios consequently requires a combination of different collider strate- 
gies in order to disentagle the larger number of possible signals. Analogously to the previous 
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Figure 10: Missing energy distribution for the events (signal+background) including 
two leptons in the final state expected for the models xen2-l,. . . ,xen2-5 sampled from 
the chains relative to benchmark xen2. This result has been obtained in a run of the 
LHC with an energy in the center of mass of 7 TeV and a luminosity of 4.7 fb _1 . The 
left panel refers to the signal region dubbed MT2 while the right panel corresponds to the 
signal region OS defined in Ref. [63]. 



benchmark we still stick on the direct production of EW interacting particles. In addition 
to three-lepton signals, however, we have also employed detection strategies of events fea- 
turing two leptons in the final state [63] which improve the sensitivity of the previous search 
in the case of the presence of sleptons in the decay chain of X\ an d are, moreover, sensitive 
to direct slepton production as well as other gaugino production processes like X\X\ which 
can be potentially relevant. 

The outcome of our analysis, performed on 5 representative points (denoted by xen2-l, 
. . . ,xen2-5) sampled from the posterior PDF of xen2 and reported in Tabs. (4) and (5), 
shows that the collider prospects are mainly influenced by the value of M2 and by the slepton 
mass scales while the [i parameter, partially fixed by DD, is important in determining the 
composition of X\ an d X*2- ^ * s a ^ so evident that the analysis employed are not sensitive 
to low mass splittings among the relevant particles. Results of a LHC run simulated from 
these points focused on final states including two and three leptons are shown in Fig. 10. 

We conclude this subsection commenting again on the possible presence of light third 
family squarks. The higher value of the DM mass (compared to xenl) weaken sensitively 
the current limits, then in principle they could account for relevant signals. However in 
the setup we are considering they do not impact the DM relic density since they cannot 
be enough light to give rise to coannihilations (concerning this claim, see also the outcome 
of the next subsection). For these reasons we will not consider this possibility since our 
principal aim is to investigate MSSM realizations which enforce correlations among the 
different experimental observables. Moreover, we notice that also in the present benchmark 
scenario we are exploring regions of the parameter space where the SI cross section is 
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Figure 11: Points of the Markov chain associated to the benchmark xen3, selected ac- 
cording the criteria depicted in the text and passing flavor and relic density constraints. 
In each panel the points are reported in the plane (m x , Clh 2 ). In the upper left panel 
the points are colored according the corresponding value of the DM higgsino fraction. 
In the upper right panel instead the color pattern is determined by the quantity Aa = 
(mA — iri x )/mA- The bottom left and bottom right panels are obtained in analogous man- 
ner but the quantities discriminated by the colors are respectively A& = (m^ — m x ) /m^ 
and At = (m^ — m x )/rni . 



dominated by the parameters of the Higgs sector. In absence of further constraints (for 
example coming from the relic density), the scan on the mass of the third family squarks 
(performed in the DD analysis) is uniquely influenced by the bound on the Higgs mass 
implemented in the Likelihood function which tends to push them towards high values (see 
also the next section for further details on this aspect). 

4.3 Analysis of xen3 

This benchmark is designed to probe coannihilation scenarios with third family squarks 
very close to the current experimental limits relative to their direct collider production. 
As already mentioned, the contribution to the SI cross section of the squark mediated 
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Table 6: Summary table of the points relative to the scenario xen3 used for the collider 
analysis. The last column reports the main decay channels of the lightest stop and sbot- 
tom in the models considered. 

interactions is enhanced for low DM-squark mass splittings and it dominates the cross- 
section itself in the limit of very low higgsino fraction. At the same time these low mass 
splittings enforce coannihilations driving the DM relic density towards its correct value and 
imply very simple and definite decay chains at collider. 

As for the other benchmark points, we start the analysis of xen3 reprocessing the 
Markov chains produced by the DD simulation applying to the sampled models the current 
bounds from the DM relic density and flavor physics. The points passing these bounds 
are shown in Fig. 11, where the color patterns follow the higgsino fraction (top left panel), 
(m-A — m x ) /uia (top right panel), (bottom left panel) and, finally A( (bottom right 
panel). The outcome of this analysis shows that extremely pure higgsino DM configu- 
rations are favored by the interplay of the DD simulation and the additional constraints 
subsequently imposed on the resulting Markov chains. Moreover, bino-like DM configu- 
rations passing these additional constraints (and corresponding to the original benchmark 
scenario) achieve the correct relic density mainly through the Higgs funnel mechanism while 
coannihilations are disfavored by the requirements on the Higgs mass (a sizable fraction of 
points in the chains is also excluded by the bounds from 6—7-57 and B s — > This 
outcome confirms the difficulty already mentioned in section 2, of realizing configurations 
in which the SI cross section is dominated by squark-mediated interactions and thus points 
again towards the need of combining the information coming from DD with the fit of some 
other observable, like for instance the relic density, which selects more efficiently the coan- 
nihilation scenario. 

The different scenarios emerging from Fig. 11 and corresponding to bino-like DM and 
higgsino-like DM candidates thermally produced in various ways (not only through coan- 
nihilations with light squarks) have different collider signatures. Since therefore in light of 
our DD simulation we can not associate to the present benchmark scenario a single clean 
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Table 7: Values of the non-SM visible cross sections of the points used for the collider 
analysis with respect to the 4 signal regions defined in Ref. [67]. This analysis is aimed 
at exploring processes leading to sensitively soft jets and then is particularly suitable to 
probe coannihilation scenarios. 

prediction for the LHC, as for instance in the case of xenl, we will consider in the following 
different cases separately, arguing that, in spite of the already mentioned difficulties, the 
coannihilation setup involving a binodike DM candidate and a stop or a sbottom NLSP 
has different collider prospects compared to, for instance, configurations corresponding to 
higgsinodike DM candidates. Analogously to the case of xen2, we will restrict to regions 
of the (M x ,asi) plane exhibiting a high statistical weight and capturing all the features 
relevant for the present analysis, namely a region corresponding to DM masses lower than 
250 GeV. 

We now focus on the binodike configurations emerging from our Monte Carlo scans. 
Given the DM mass range relevant for the present analysis, production processes involving 
gluinos and first generation squarks are still severely constrained. If stops and sbottoms 
lie in the TeV range, like in the CP-odd resonance scenario, the dominant process in this 
case is again gaugino production with typical cross section of only few fb and thus very 
difficult to detect in the next future. For this reason we will not further consider this 
possibility in what follows. On the contrary the coannihilation scenario involving sbottom 
and squarks NLSP features sizable production rates for these particles. A coannihilating 
sbottom NLSP can only decay through the channel b — > bx± while in the case of a stop 
NLSP the dominant decay channel is often t — > ex?, arising at onedoop, while the other 
possible decay is £ — 5- bff'xi, with / and /' being two light fermions. These decays, 
however, evade most of the current SUSY searches since their products are too soft to 
pass the analysis cuts on the px, or even the requirements for correct reconstruction of the 
final states. Such configurations thus require search strategies optimized for very low mass 
splittings or alternatively can be probed using initial state radiation. In this work we have 
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Table 8: Values of the non-SM visible cross sections of the points used for the collider 
analysis with respect to the 3 signal regions at 2-jet minimal multiplicity defined in 
Ref. [69]. This analysis is aimed at exploring decay chains with low mass splittings. 

focused on the case of the initial radiation of a jet employing the dedicated ATLAS search 
of Ref. [81]. Good prospects of probing scenarios with very low mass splitting have been 
obtained as well from the study of monophoton events in Ref. [74]. 

To this purpose we have sampled from the PDF obtained from our DD analysis a 
few representative models (listed in Tab. 6) and generated from them LHC final states 
involving pairs of stop and sbottom squark with an additional jet. The results, reported 
in Tab. 7, show that this coannihilation scenario is rather efficiently probed, with some 
points already excluded, especially at low masses while we have a significant decrease in the 
performance at higher values since the stop/sbottom production cross sections are rather 
steeply decreasing functions of the masses of these particles. These results further stress 
the importance of increasing the performance in reconstructing the DM mass (in this case 
related by coannihilation requirements to the squark ones) since even slightly different 
determinations drastically alter the collider detection prospects. 

The results of monojet searches can be interpreted in a large variety of models including 
some non-SUSY extension of the SM, as for instance the one considered in Ref. [81]. In 
addition to the monojet analysis, we have then employed a search based on two jets in 
the final state (see also Ref. [82]) featuring a signal region referred as A' [69] specifically 
designed for being sensitive to low particle mass splittings 11 . As evident from Tab. 8, this 
kind of search has shown a good sensitivity for some of the sample points although it turned 
out to be in general less effective than the monojet analysis. 

We conclude this section focusing on the collider signatures of the higgsino DM scenarios 

A preliminar analysis updated to the 2012 data set with 8 TeV of center of mass energy is already 
available [70]. However a slightly different set of cuts have been employed. The result here shown are 
related to the 7 TeV data set, no sensitive improvement have been found once performing the update in 
luminosity and center of mass energy. 
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Table 9: Mass of the lightest neutralino and chargino for the additional points at hig- 
gsino DM for which no additional requirment on the rest of the spectrum is imposed. 
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Table 10: Values of the non-SM visible cross section of the points reporterd in tab. (9) 
relative to the 4 signal regions designed for the collider search described in [67]. 

emerging from our Monte Carlo scans. The main peculiar feature of these SUSY spectra is 
the presence of a chargino and a second neutralino almost degenerate with the LSP, with 
the masses of all the three particles controlled by the [i parameter. The correct DM relic 
density is guaranteed by the higgsino nature of the DM candidate, so no other requirement 
is imposed on the rest of the spectrum. If this is the case, collider phenomenology is 
characterized again by neutralino/chargino production (namely the production of %\ X% but 
also XiXi an(I XiXi) as follows by the mass degeneracy of the three lightest higgsinos, thus 
featuring a possible signal into monojets. We have accordingly simulated the production 
of a pair of lightest neutralino/charginos with one additional jets for a sample drawn from 
our chains of eight points (see Tab. 9) with mass as low as 134 GeV in order to take 
into account the outcome of the previous section as well. The values of the visible cross 
sections are however significantly lower (see Tab. 10) with respect to the previously discussed 
coannihilation scenario which is then better probed by mono jet searches. 
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Figure 12: Points in the Markov chain associated to the benchmark xen4 passing fla- 
vor and relic density constraints. The points in the two plots follow color patterns de- 
termined, respectively, by higgsino fraction and A^. By comparing the two panels one 
notices that a large fraction of points at gaugino DM clusters around the value ~ 0.5. 

4.4 Analysis of xen4 

The last benchmark point considered in the analysis has been also designed, similarly to 
xen3, to investigate a scenario where DM is produced through coanniliations involving 
squarks of the third family. In this case, however, the DM mass is significantly higher 
than in the case of the previous benchmark. Fig. 12 clearly shows that no coannihilation 
configurations emerged from the interplay of the DD simulations and the application of the 
relic density and flavor bounds, being the PDF samples strongly dominated by higgsino 
DM candidates and, to some extent, by CP-odd Higgs funnel configurations. In view of 
this outcome, we regard the reconstruction of this point as unsuccessful and thus we will 
not pursue any collider analysis for this scenario. 

5 Conclusions 

In this paper we have explored a new strategy to combine direct and collider dark matter 
searches: employing the potential of the upcoming generation of 1-ton direct detection 
experiments, we have studied the possibility of translating the information encoded in an 
hypothetically discovered direct detection signal into certain classes of expected signals at 
the LHC. 

The analysis has been performed in the context of the so-called phenomenological Min- 
imal Supersymmetric Standard Model where a large number of possible supersymmetric 
configurations is conveniently described by a limited number of parameters. We focused 
on four benchmark points representative of WIMP dark matter candidates which should 
be soon probed by the upcoming generation of 1-ton direct detection experiments and 
corresponding to different dark matter production scenarios. In all cases a clear correla- 
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tion pattern between direct detection and LHC observables is enforced by the underlying 
mechanism accounting for the dark matter relic density 

For each benchmark, we tested the capability of a 1-ton direct detection experiment 
resembling the XEN0N1T experiment in reconstructing the dark matter mass and spin- 
independent cross section within the framework of a Bayesian analysis supported by a 
Markov Chain Monte Carlo scanning technique. In our analysis we carefully addressed the 
impact on the reconstructions of the astrophysical uncertainties entering the computation 
of the differential count rates, of possible statistical limitations related to the finite number 
of measured recoil energies as well as of the choice of the priors for the model parameters. 

We derived through this procedure the accuracy within which the dark matter mass 
and scattering cross section are reconstructed by the Xenon experiment in its 1-ton configu- 
ration. We then focused on the regions of the plane dark matter mass versus scattering cross 
section enclosing the 95% of the posterior probability. Within these regions we selected a 
subset of models by reprocessing our Markov chains introducing additional constraints from 
the relic density, the anomalous magnetic moment of the muon and B-physics processes. 
For some of these models we simulated the ATLAS detector response in its configuration 
with 7 TeV center of mass energy and carefully studied the prospects of detecting certain 
classes of final states for which dedicated search strategies have been developed by the 
ATLAS collaboration. 

The procedure proposed here defines a method to efficiently exploit the upcoming gen- 
erations of 1-ton direct detection experiments to forecast valuable predictions for the LHC 
and it constitutes the main achievement of this paper. The results obtained within this 
framework are particularly promising for the benchmark xenl which is representative of a 
light bino-like neutralino produced via resonant annihilations. In this case a definite pre- 
diction for a LHC final state including three leptons and missing energy can be elaborated 
starting from the informations derived from the simulated 1-ton direct detection experi- 
ment. The other benchmarks, instead, are more sensitive to the direct detection statistical 
limitations discussed in the text and affecting the accuracy of the dark matter mass recon- 
struction. Larger exposures in the direct detection (here we considered 1000 kg-years) and 
eventually more informative Likelihoods can however further enhance the effectiveness of 
our approach and allow to draw interesting conclusions also in the case of challenging dark 
matter candidates as for instance xen4. 

In summary, we have proposed a strategy to systematically translate an hypothetical 
direct detection signal into a prediction for the LHC. This method has the potential to 
significantly strengthen the complementarity between the direct detection technique and 
dark matter searches at collider and will be further investigated in future works. 

Acknowledgments 

It is a pleasure to thank Antonio Masiero for his collaboration during the early stages of the 
project and useful discussions on the main idea behind this work. We also thank Umberto 
de Sanctis, Alberto Tonero, Michele Pinamonti, Andrea De Simone and Florian Bonnet for 
fruitful conversations. These results have been partially obtained making use of the bwGRiD 



- 35 - 



Cluster (http://www.bw-grid.de), member of the German D-Grid initiative, funded by the 
Ministry for Education and Research (Bundesministerium fur Bildung und Forschung) and 
the Ministry for Science, Research and Arts Baden- Wiirttemberg (Ministerium fur Wis- 
senschaft, Forschung und Kunst Baden- Wiirttemberg). G.A and P.U. acknowledges partial 
support from the European Union FP7 ITN-INVISIBLES (Marie Curie Actions, PITN- 
GA-2011- 289442). G.A also thanks the Galileo Galilei Institute for the warm hospitality 
during part of the realization of this work. 

References 

[1] P. J. Fox, R. Harnik, J. Kopp, and Y. Tsai, "Missing Energy Signatures of Dark Matter at 
the LHC," Phys.Rev., vol. D85, p. 056011, 2012, 1109.4398. 

[2] et ai, "Dark Matter Results from 225 Live Days of XENON100 Data," 2012, 1207.5988. 

[3] R. Catena and P. Ullio, "A novel determination of the local dark matter density," JCAP, 
vol. 1008, p. 004, 2010, 0907.0018. 

[4] R. Catena and P. Ullio, "The local dark matter phase-space density and impact on WIMP 
direct detection," JCAP, vol. 1205, p. 005, 2012, 1111.3556. 

[5] F. Iocco, M. Pato, G. Bertone, and P. Jetzer, "Dark Matter distribution in the Milky Way: 
microlensing and dynamical constraints," JCAP, vol. 1111, p. 029, 2011, 1107.5810. 

[6] S. Garbari, J. I. Read, and G. Lake, "Limits on the local dark matter density," 
Mon.Not.Roy.Astron.Soc., vol. 416, pp. 2318-2340, 2011, 1105.6339. 

[7] M. Pato, O. Agertz, G. Bertone, B. Moore, and R. Teyssier, "Systematic uncertainties in the 
determination of the local dark matter density," Phys.Rev.. vol. D82, p. 023531, 2010, 
1006.1322. 

[8] P. Salucci, F. Nesti, G. Gentile, and C. Martins, "The dark matter density at the Sun's 
location," Astron.Astrophys., vol. 523, p. A83, 2010, 1003.3101. 

[9] W. de Boer and M. Weber, "The Dark Matter Density in the Solar Neighborhood 
reconsidered," JCAP, vol. 1104, p. 002, 2011, 1011.6323. 

[10] S. H. Hansen, B. Moore, M. Zemp, and J. Stadel, "A Universal velocity distribution of 
relaxed collisionless structures," JCAP, vol. 0601, p. 014, 2006, astro-ph/0505420. 

[11] C. Strege, R. Trotta, G. Bertone, A. H. Peter, and P. Scott, "Fundamental statistical 
limitations of future dark matter direct detection experiments," Phys.Rev., vol. D86, 
p. 023507, 2012, 1201.3631. 

[12] A. Djouadi et at, "The Minimal supersymmetric standard model: Group summary report," 
1998, hep-ph/9901246. 

[13] G. Aad et ai, "Search for squarks and gluinos using final states with jets and missing 

transverse momentum with the ATLAS detector in sqrt(s) = 7 TeV proton-proton collisions," 
2011, 1109.6572. 

[14] S. Chatrchyan et ai, "Search for Supersymmetry at the LHC in Events with Jets and 
Missing Transverse Energy," Phys.Rev. Lett., vol. 107, p. 221804, 2011, 1109.2352. 

[15] A. Djouadi, J.-L. Kneur, and G. Moultaka, "SuSpect: A Fortran code for the 



- 36 - 



supersymmetric and Higgs particle spectrum in the MSSM," Comput.Phys.Commun., 
vol. 176, pp. 426-455, 2007, hep-ph/0211331. 

[16] G. Aad et ai, "Observation of a new particle in the search for the Standard Model Higgs 
boson with the ATLAS detector at the LHC," Phys.Lett., vol. B716, pp. 1-29, 2012, 
1207.7214. 

[17] P. P. Giardino, K. Kamiike, M. Raidal, and A. Strumia, "Is the resonance at 125 GeV the 
Higgs boson?," 2012, 1207.1347. 

[18] D. Carmi, A. Falkowski, E. Kufiik, T. Volansky, and J. Zupan, "Higgs After the Discovery: A 
Status Report," 2012, 1207.1718. 

[19] H. E. Haber, "Challenges for nonminimal Higgs searches at future colliders," 1995, 
hep-ph/9505240. 

[20] H. Baer, V. Barger, and A. Mustafayev, "Implications of a 125 GeV Higgs scalar for LHC 
SUSY and neutralino dark matter searches," Phys.Rev., vol. D85, p. 075010, 2012, 1112.3017. 

[21] A. Arbey, M. Battaglia, A. Djouadi, F. Mahmoudi, and J. Quevillon, "Implications of a 125 
GeV Higgs for supersymmetric models," Phys.Lett., vol. B708, pp. 162-169, 2012, 1112.3028. 

[22] L. J. Hall, D. Pinner, and J. T. Ruderman, "A Natural SUSY Higgs Near 126 GeV," JEEP, 
vol. 1204, p. 131, 2012, 1112.2703. 

[23] A. Arbey, M. Battaglia, A. Djouadi, and F. Mahmoudi, "The Higgs sector of the 

phenomenological MSSM in the light of the Higgs boson discovery," 2012, 1207.1348. 

[24] G. Aad et ai, "Search for neutral MSSM Higgs bosons decaying to t + t~ pairs in 

proton-proton collisions at sqrt(s) = 7 TeV with the ATLAS detector," Phys.Lett., vol. B705, 
pp. 174-192, 2011, 1107.5003. 

[25] S. Chatrchyan et al, "Search for neutral Higgs bosons decaying to tau pairs in pp collisions 
at sqrt(s)=7 TeV," Phys.Lett., vol. B713, pp. 68-90, 2012, 1202.4083. 

[26] P. P. Giardino, K. Kamiike, M. Raidal, and A. Strumia, "Reconstructing Higgs boson 
properties from the LHC and Tevatron data," JEEP, vol. 1206, p. 117, 2012, 1203.4254. 

[27] T. Moroi, "The Muon anomalous magnetic dipole moment in the minimal supersymmetric 
standard model," Phys.Rev., vol. D53, pp. 6565-6575, 1996, hep-ph/9512396. 

[28] L. Calibbi, R. Hodgkinson, J. Jones Perez, A. Masiero, and O. Vives, "Flavour and Collider 
Interplay for SUSY at LHC7," Eur.Phys.J., vol. C72, p. 1863, 2012, 1111.0176. 

[29] E. Lunghi and J. Matias, "Huge right-handed current effects in B — > K * (Ktt)1 + E in 
supersymmetry," JEEP, vol. 0704, p. 058, 2007, hep-ph/0612166. 

[30] F. Mahmoudi, "Superlso: A Program for calculating the isospin asymmetry of B — > K * 7 in 
the MSSM," Comput.Phys.Commun., vol. 178, pp. 745-754, 2008, 0710.2067. 

[31] F. Mahmoudi, "Superlso v2.3: A Program for calculating flavor physics observables in 
Supersymmetry," Comput.Phys.Commun., vol. 180, pp. 1579-1613, 2009, 0808.3144. 

[32] S. Chatrchyan et ai, "Search for B° s tofi + [T and 5°fo^+/i~ decays," JEEP, vol. 1204, p. 033, 
2012, 1203.3976. 

[33] R. Aaij et at, "Strong constraints on the rare decays B s — > and B — > 

Phys.Rev.Lett., vol. 108, p. 231801, 2012, 1203.4493. 

[34] P. Gondolo et at, "DarkSUSY: Computing supersymmetric dark matter properties 



-37- 



numerically," J GAP, vol. 0407, p. 008, 2004, astro-ph/0406204. 

[35] A. Hryczuk, R. Iengo, and P. Ullio, "Relic densities including Sommerfeld enhancements in 
the MSSM," JEEP, vol. 1103, p. 069, 2011, 1010.2172. 

[36] N. Arkani-Hamed, A. Delgado, and G. Giudice, "The Well-tempered neutralino," Nucl.Phys., 
vol. B741, pp. 108-130, 2006, hep-ph/0601041. 

[37] J. R. Ellis, K. A. Olive, Y. Santoso, and V. C. Spanos, "Update on the direct detection of 
supersymmetric dark matter," Phys.Rev., vol. D71, p. 095007, 2005, hep-ph/0502001. 

[38] M. Drees and M. M. Nojiri, "New contributions to coherent neutralino - nucleus scattering," 
Phys.Rev., vol. D47, pp. 4226-4232, 1993, hep-ph/9210272. 

[39] M. Drees and M. Nojiri, "Neutralino - nucleoli scattering revisited," Phys.Rev., vol. D48, 
pp. 3483-3501, 1993, hep-ph/9307208. 

[40] J. Binney and S. Tremaine, "Galactic Dynamics: Second Edition," Princeton University 
Press, 2008. 

[41] H. Freudenreich, "Cobe's galactic bar and disk," Astrophys.J., vol. 492, pp. 495-510, 1998, 
astro-ph/9707340. 

[42] T. M. Dame, "The distribution of neutral gas in the milky way," AIP Conference 
Proceedings, vol. 278, no. 1, pp. 267-278, 1992. 

[43] J. F. Navarro, E. Hayashi, C. Power, A. Jenkins, C. S. Frenk, et ai, "The Inner structure of 
Lambda-CDM halos 3: Universality and asymptotic slopes," Mon.Not.Roy.Astron.Soc, 
vol. 349, p. 1039, 2004, astro-ph/0311231. 

[44] A. W. Graham, D. Merritt, B. Moore, J. Diemand, and B. Terzic, "Empirical Models for 
Dark Matter Halos. III. The Kormendy relation and the log(rho e ) — log(R e ) relation," 
Astron.J., vol. 132, pp. 2711-2716, 2006, astro-ph/0608614. 

[45] G. Bryan and M. Norman, "Statistical properties of x-ray clusters: Analytic and numerical 
comparisons," Astrophys.J., vol. 495, p. 80, 1998, astro-ph/9710107. 

[46] E. Komatsu et ai, "Seven- Year Wilkinson Microwave Anisotropy Probe (WMAP) 
Observations: Cosmological Interpretation," Astrophys.J. Supply vol. 192, p. 18, 2011, 
1001.4538. 

[47] X. Xue et ai, "The Milky Way's Circular Velocity Curve to 60 kpc and an Estimate of the 
Dark Matter Halo Mass from Kinematics of 2400 SDSS Blue Horizontal Branch Stars," 
Astrophys.J., vol. 684, pp. 1143-1158, 2008, 0801.1232. 

[48] Y. Akrami, C. Savage, P. Scott, J. Conrad, and J. Edsjo, "How well will ton-scale dark 

matter direct detection experiments constrain minimal supersymmetry?," JCAP, vol. 1104, 
p. 012, 2011, 1011.4318. 

[49] G. Bertone, D. G. Cerdeno, M. Fornasa, R. Ruiz de Austri, C. Strege, et ai, "Global fits of 
the cMSSM including the first LHC and XENON100 data," JCAP, vol. 1201, p. 015, 2012, 
1107.1715. 

[50] D. G. Cerdeno, M. Fornasa, J.-H. Huh, and M. Peiro, "Nuclear uncertainties in the 
spin-dependent structure functions for direct dark matter detection," 2012, 1208.6426. 

[51] Y. Akrami, C. Savage, P. Scott, J. Conrad, and J. Edsjo, "Statistical coverage for 

supersymmetric parameter estimation: a case study with direct detection of dark matter," 
JCAP, vol. 1107, p. 002, 2011, 1011.4297. 



- 38 - 



[52] R. Trotta, F. Feroz, M. P. Hobson, L. Roszkowski, and R. Ruiz de Austri, "The Impact of 
priors and observables on parameter inferences in the Constrained MSSM," JHEP, vol. 0812, 
p. 024, 2008, 0809.3792. 

[53] G. Bertone, D. G. Cerdeno, M. Fornasa, R. R. de Austri, and R. Trotta, "Identification of 
Dark Matter particles with LHC and direct detection data," Phys.Rev., vol. D82, p. 055008, 
2010, 1005.4280. 

[54] G. B. Gehnini, P. Gondolo, A. Soldatenko, and C. E. Yaguna, "Direct detection of neutralino 
dark mattter in non- standard cosmologies," Phys. Rev., vol. D76, p. 015010, 2007, 
hep-ph/0610379. 

[55] B. S. Acharya, G. Kane, S. Watson, and P. Kumar, "A Non-thermal WIMP Miracle," Phys. 
Rev., vol. D80, p. 083529, 2009, 0908.2430. 

[56] G. Arcadi and P. Ullio, "Accurate estimate of the relic density and the kinetic decoupling in 
non-thermal dark matter models," Phys.Rev., vol. D84, p. 043520, 2011, 1104.3591. 

[57] R. Catena, N. Fornengo, A. Masiero, M. Pietroni, and F. Rosati, "Dark matter relic 
abundance and scalar - tensor dark energy," Phys.Rev., vol. D70, p. 063519, 2004, 
astro-ph/0403614. 

[58] R. Catena, N. Fornengo, A. Masiero, M. Pietroni, and M. Schelke, "Enlarging mSUGRA 
parameter space by decreasing pre-BBN Hubble rate in Scalar- Tensor Cosmologies," JHEP, 
vol. 0810, p. 003, 2008, 0712.3173. 

[59] G. Aad et ai, "Expected Performance of the ATLAS Experiment - Detector, Trigger and 
Physics," 2009, 0901.0512. 

[60] J. Alwall, M. Herquet, F. Maltoni, O. Mattelaer, and T. Stelzer, "MadGraph 5 : Going 
Beyond," JHEP, vol. 1106, p. 128, 2011, 1106.0522. 

[61] W. Beenakker, R. Hopker, and M. Spira, "PROSPINO: A Program for the production of 
supersymmetric particles in next-to-leading order QCD," 1996, hep-ph/9611232. 

[62] G. Aad et al., "Search for direct production of charginos and neutralinos in events with three 
leptons and missing transverse momentum in sqrt(s) = 7 TeV pp collisions with the ATLAS 
detector," 2012, 1208.3144. 

[63] G. Aad et al., "Search for direct slepton and gaugino production in final states with two 
leptons and missing transverse momentum with the ATLAS detector in pp collisions at 
sqrt(s) = 7 TeV," 2012, 1208.2884. 

[64] G. Aad et al., "Search for top and bottom squarks from gluino pair production in final states 
with missing transverse energy and at least three b-jets with the ATLAS detector," 2012, 
1207.4686. 

[65] G. Aad et al., "Search for direct top squark pair production in final states with one isolated 
lepton, jets, and missing transverse momentum in sqrt(s) = 7 TeV pp collisions using 4.7 fb-1 
of ATLAS data," 2012, 1208.2590. 

[66] G. Aad et al., "Search for a supersymmetric partner to the top quark in final states with jets 
and missing transverse momentum at sqrt(s) = 7 TeV with the ATLAS detector," 2012, 
1208.1447. 

[67] "Search for dark matter candidates and large extra dimensions in events with a jet and 

missing transverse momentum with the atlas detector," Tech. Rep. ATLAS-CONF-2012-084, 



- 39 - 



CERN, Geneva, Jul 2012. 

[68] "Search for scalar bottom pair production in E>nal states with missing transverse momentum 
and two b-jets in pp collisions at As=7 tev with the atlas detector," Tech. Rep. 
ATLAS-CONF-2012-106, CERN, Geneva, Aug 2012. 

[69] G. Aad et al., "Search for squarks and gluinos with the ATLAS detector in final states with 
jets and missing transverse momentum using 4.7 fb~ 1 of sqrt(s) = 7 TeV proton-proton 
collision data," 2012, 1208.0949. 

[70] "Search for squarks and gluinos with the atlas detector using final states with jets and 

missing transverse momentum and 5.8 fb _1 of y/s=8 tev proton-proton collision data," Tech. 
Rep. ATLAS-CONF-2012-109, CERN, Geneva, Aug 2012. 

[71] H. Baer, V. Barger, A. Lessa, W. Sreethawong, and X. Tata, "Wh plus missing- S^signature 
from gaugino pair production at the LHC," Phys.Rev., vol. D85, p. 055022, 2012, 1201.2949. 

[72] J. A. Conley, J. S. Gainer, J. L. Hewett, M. P. Le, and T. G. Rizzo, "Supersymmetry 
Without Prejudice at the 7 TeV LHC," Physical Review D, 2011, 1103.1697. 

[73] H. K. Dreiner, M. Kramer, and J. Tattersall, "How low can SUSY go? Matching, monojets 
and compressed spectra," Europhys.Lett., vol. 99, p. 61001, 2012, 1207.1613. 

[74] G. Belanger, M. Heikinheimo, and V. Sanz, "Model-Independent Bounds on Squarks from 
Monophoton Searches," JHEP, vol. 1208, p. 151, 2012, 1205.1463. 

[75] G. Aad et al., "Search for scalar bottom pair production with the ATLAS detector in pp 
Collisions at sqrts = 7 TeV," Phys.Rev. Lett., vol. 108, p. 181802, 2012, 1112.3832. 

[76] G. Aad et al., "Search for light top squark pair production in final states with leptons and 
b-jets with the ATLAS detector in sqrt(s) = 7 TeV proton-proton collisions," 2012, 1209.2102. 

[77] A. Choudhury and A. Datta, "Many faces of low mass neutralino dark matter in the 

unconstrained MSSM, LHC data and new signals," JHEP, vol. 1206, p. 006, 2012, 1203.4106. 

[78] O. Adriani et al., "PAMELA results on the cosmic-ray antiproton flux from 60 MeV to 180 
GeV in kinetic energy," Phys.Rev. Lett., vol. 105, p. 121101, 2010, 1007.0821. 

[79] C. Weniger, "A Tentative Gamma-Ray Line from Dark Matter Annihilation at the Fermi 
Large Area Telescope," JCAP, vol. 1208, p. 007, 2012, 1204.2797. 

[80] W. Buchmuller and M. Garny, "Decaying vs Annihilating Dark Matter in Light of a 
Tentative Gamma-Ray Line," JCAP, vol. 1208, p. 035, 2012, 1206.7056. 

[81] G. Aad et al., "Search for dark matter candidates and large extra dimensions in events with 
a jet and missing transverse momentum with the ATLAS detector," 2012, 1210.4491. 

[82] A. Choudhury and A. Datta, "New limits on top squark NLSP from LHC 4.7 fb^ 1 data," 
2012, 1207.1846. 



-40- 



